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Graphene has been considered to be a promisine; material for future electron- 

n n 

ics due to its extraordinary properties such as ^igh carrier mobility [1, 2j, ther- 



mal conductivity [3], and strong break strength j^. Although the extremely high 
electrical conductivity makes graphene a candidate for replacing silicon-based 
electronics, Klein tunneling causes that electrical transport of Dirac fermions 
is insensitive to electrostatic potentials, resulting a low current on/off ratio 
of graphene-based field effect transistors |5^-Q]. In order to increase the cur- 
rent on/off" ratio, vertical graphene heterostructures have been introduced as 
an alternative device architecture!^, Q' by using quantum tunneling. Here, 
in addition to the experimental observations which gave us a large current 
on/off ratio of graphene/MoS2 heterostructures, we present that the current 
on/off" ratio is enhanced by as large as two orders of magnitude compared 
to the previous report [sj]. Moreover, we also report a novel utility of the 
graphene/MoS2 heterostructures as a spin-filter. When a thin M0S2 layer sand- 
wiched between two graphene sheets becomes magnetic because of the inversion 
symmetry breakingjiot], Dirac fermions with diff"erent spins feel different heights 
of a tunnel barrier, resulting spin-dependent tunneling. Our findings will de- 
velop the present graphene heterostructures for electronic devices by improving 
the device performance and by adding the functionality to the existing graphene 
heterostructures. 

Very recently, increasing interest has been focussed on an alternative graphene device 
structure by using quantum tunneling. In a graphene/silicon junction, large current on/off 
ratio was achieved by controlling the Schottky barrier at the interface However, for 
graphene deposited on silicon substrate, carrier mobility is expected to decrease because of 



carrier inhomogeneity induced by the substrate [12j. Meanwhile, L. Britnell et al. reported 
the possibility of a graphene field-effect-transistor based on vertical heterostructures with 
atomically thin insulating barriers, hexagonal boron nitride (/iBN) and molybdenum disul- 



fide (MoS2)[8|. /iBN has gained burgeoning interest as a material for use in graphene devices 
because the encapsulation of graphene by /iBN maintains the high electronic quality of pris- 



tine one|l3l-ll5|. In spite of this advantage, the large band gap of hBN (~ 5.97 e\^[16i]) causes 
an insufficient current on/off ratio. Instead of /iBN, the larger on/off ratio was observed in 
graphene heterostructure with M0S2, owing to its smaller bandgap than hBN. Therefore, the 



graphene heterostructure with M0S2 can be an important building block of graphene-based 
electronics, and it is natural to investigate possible functional devices utilizing its advantages 
for applications. 

Herein, we present not only the improvement of the current on/off ratio of the existing 



graphene/MoS2 heterostructuresjs] but also an application of the heterostructure in spin- 
tronics by producing spin-dependent transport. First, we show that there emerge current 
peaks for a graphene/MoS2/grarphene nanoribbon (GNR) heterostructure. This finding has 
potential for the use of the current peaks, resulting in the improvement of the current on/off 



ratio. Second, the existence of magnetic properties in few-layer MoS9|17l-l22| can lead to 
the spin-polarized current in the graphene heterostructures. We show that the graphene 
heterostructure can be a perfect spin-filter for holes with the electron-hole asymmetric spin 
splitting of M0S2 lo|. Such tunneling phenomena may offer further potential applications 
in graphene-based electronics and spintronics. 

Now, we consider a heterostructure, which consists of an atomically thin M0S2 layer 
sandwiched by two graphene sheets as shown in Fig. [1^. It is well known that few-layer 
M0S2 is an insulator with finite bandgaps; ~ 1.9 eV direct bandgap near the K- valley and 



~ 1.2 — 1.4 eV indirect bandgap depending on the number of M0S2 layers 23|]. The M0S2 
layer of the heterostructure becomes a tunnel barrier for Dirac fermions, and both graphene 
sheets play the role of high-quality source and drain electrodes. Dirac fermions experience 
the direct bandgap near K- valley of M0S2 because of the momentum conservation, neglecting 
electron-phonon scattering processes. 

We can obtain the tunneling current through the M0S2 insulating barrier as below 

/+00 
Ds {E, Vh) Dd {E, H) T {E) [fs {E, H, Vg) - fn {E, Vg)] dE, (1) 
-00 

where jo = (qvp) / {2TrLl) is the unit of current density with electric charge of carriers, q, and 
the characteristic length of the system, Lq. We can calculate the transmission probability, 
T{E), with a quantum- mechanical consideration (see Supplementary Information). Here, 
Di and /j are density of states of graphene and the Fermi-Dirac distribution where i = 
S, D represent source and drain graphene electrodes on both sides of the M0S2 layer, 
respectively. We assume that the gate voltage induces the equal carrier concentration in 
both graphene sheets; n = aVg where a = 6.16 x 10^^ Cs'^kg~'^mr^ in a case of a 300 nm 
thick Si02 substrate. For a given Vg, chemical potentials in both graphene sheets are equally 




FIG. 1. Graphene heterostructures with M0S2 thin layers, a, b, c, Each panel shows ener- 
getic diagrams for quantum tunnehng through M0S2 insulating barriers for different heterostruc- 
tures; graphene/MoS2 /graphene, graphene/MoS2/AGNR, and graphene/m-MoS2 /graphene. Since 
the charge neutral point of graphene is asymmetrically laid, electrons and holes experience different 
tunnel barriers, Ac and A^,, respectively, b, If one graphene electrode is replaced by a narrow GNR, 
density of states is changed, showing the one-dimensional nature, c. When the M0S2 is magnetic, 
the barrier height becomes spin-dependent, with different spin-splitting energies 6c and 5^. d, e, 
f. Schematic diagrams of different heterostructures corresponding to a, b, c, respectively. M0S2 is 
sandwiched by two graphenes which act as source and drain electrodes. The equilibrium chemical 
potential, /Uq, is induced via back gate, and the tunneling current is generated by source-drain bias 
voltage, Vb- 

given as yU = Hvp^Jt^ \n\ = hvF\/T^a \Vg\. The equal chemical potentials produce no net 
tunneling current through the M0S2 barrier. Applying the bias voltage, Vf,, between the 
source and drain graphene electrodes, one can measure non-zero tunneling current through 
the heterostructure. 

Figure |2] displays the characteristics of the graphene/MoS2/graphene heterostructure. 
The calculated tunneling currents as functions of H for different Vg are shown in Fig. [2^. 
The tunneling current exhibits an exponentially increasing behaviour and becomes larger as 
Vg increases because more carriers contribute to the quantum tunneling for high Vg. The 
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FIG. 2. I-V characteristics of graphene/MoS2/graphene heterostructure. a, Tunneling 
current densities as functions of bias voltage for different gate voltages at 300 K. b, Temperature 
dependence of the tunneling current densities. Tunneling current densities are plotted in loga- 
rithmic scale as functions of gate voltage for the given bias voltage, Vb = 0.1 V. c, Asymmetric 
tunneling current density curve versus gate voltage at 300 K for the given bias voltage, = 0-1 V. 
(Inset) Energetic diagrams for electron and hole tunnelings. Electrons and holes are induced by 
applying the gate voltage, Vg. d, Log-linear curves of the electron and hole tunneling current 
densities as a function of gate voltage at 300 K for the given bias voltage, Vb = 0.1 V. In spite 
of different magnitudes of the current densities between electrons and holes, we can obtain similar 
current on/off ratios about 10^. A parameter jo = (qvp) / (27rLg) in use is the unit of current 
density. 



tunneling current versus gate voltage curves are plotted in Fig. ^jp. We find the ratio of the 
tunneling current between an off-state {Vg = V^) and an on-state (Vg = 150 V) at different 
temperatures. With increasing temperature, the magnitude of tunneling current correspond- 



ingly increases because more carriers are thermally activated, but the current on/off ratio 
decreases. While the current on/off ratio reaches up to 10^^ at 100 it goes down to 
10^ at 300 K. Despite of this decrease, the room-temperature current on/off ratio is still 
high as experimentally expected in Ref. 18. This means that the graphene/MoS2/graphene 
heterostructure can be a platform for graphene-based field-effect transistors. 

For graphene/MoS2 hybrid systems, the charge neutral point of Dirac cone is asymmetri- 



cally arranged between the conduction and valence bands near K- valley of MoS2[2l|]. Thus, 
electrons experience a lower tunnel barrier (A^ ~ 0.5 eV) than holes (A^, ~ 1.4 eV), where 
Ac and A„ are the barrier heights for electrons and holes, respectively. If the tunneling 
electrons and holes have the same energy, the transmission probability through the M0S2 
insulating barrier for electrons is greater than for holes because electrons experience the 
smaller tunnel barrier than holes. Figure |2t exhibits the difference between electron and 
hole tunneling currents as functions of gate voltage. We found that the tunneling current is 
indeed asymmetric with respect to the gate voltage polarity. However, as shown in Fig. [SJi, 
the on/off ratio of electron and hole tunneling currents is about the same although the hole 
tunneling current is suppressed. 

Next, we can enhance the current on/off ratio of the graphene/MoS2/graphene het- 
erostructure by considering size effects of a graphene electrode. In order to study the size 
effects, we replace a graphene electrode on one side of the sandwiched M0S2 layer with a 
GNR as shown in Fig. [T]d. Because the tunneling current density depends on the product 
of densities of states of source and drain graphene electrodes, (Eq. ([1])), one can expect the 
changes in the tunneling current. In this letter, we take into account a graphene nanoribbon 
with armchair edges (AGNR) of which the density of states exhibits the one-dimensional 
nature (see Supplementary Information). 

The resulting tunneling currents through a graphene/MoS2/AGNR heterostructure are 
plotted in Fig. |3^ as functions of H for different Vg. There emerge tunneling current 
peaks at specific H because of the one-dimensional nature of the AGNR. This feature of 
the graphene/MoS2/AGNR heterostructure plays a crucial role in the enhancement of the 
current on/off ratio. Figure 13b shows log-linear plots of the tunneling current densities 
as functions of Vg for different VJ,. The magnitude of the currents at a current peak (for 
Vf, ^ 0.1 V ^ blue lines) is about 10^ times larger than the background values (for VJ, ~ 0.11 V , 
red lines). If one vary VJ, from ^ 0.11 ^ to ^ 0.1 V simultaneously with adjusting V^, one 
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FIG. 3. I-V characteristics of graphene/MoS2/AGNR heterostructures. a, Tunneling 
current densities through a graphene/MoS2/AGNR heterostructure as functions of bias voltage for 
different gate voltages (arbitrary units). Some current peaks emerge at the specific bias voltages 
due to the one-dimensional nature of AGNRs. b, Comparison of the current on/off ratios in 
different cases; blue, red, and green solid lines depicted in a. Blue and red dashed curves indicate 
the current density versus the gate voltage for different fixed bias voltages, Vf, ~ 0.1 V and 0.11 V, 
respectively. Green solid curve shows the current density as a function of gate voltage, varying 
simultaneously bias voltage from 0.11 V to 0.1 V. 

can obtain an abrupt jump in the tunneling current density as large as 10^ (green lines). Due 
to the existence of the current peak, there emerges an abrupt jump of the current density as 
large as 10^, and the resulting current on/off ratio is enhanced up to 10*^. Therefore, the use 
of an AGNR instead of a graphene electrode is attractive for applications in graphene-based 
electronics, improving the performance of the graphene heterostructures with M0S2 layers. 

7 



In recent years, it was found to exhibit magnetic properties in an atomically thin M0S2 



layer due to several causes; sulfur- vacancy or zigzag-terminated grain boundaryjlOl. |17H22|. 
For instance, the broken inversion symmetry due to the sulfur-vacancy leads to a splitting 
between different spin states for few-layer M0S2, but there is no spin-splitting for bulk 
M0S2 10|. This allows us to assume that the M0S2 layer used in our heterostructures is thin 



enough to have a non-zero spin-splitting energy. In this case, the thin M0S2 layer can be 
treated as a magnetic insulator with spin-dependent barrier heights, U (z) = A + cr6c,v + 
{qVbz/d) where cr = ± represents different spins and 6c,v indicate spin-splitting energies in 
the conduction and valence bands of M0S2, respectively. We focus on the spin-splitting 
near K- valley of M0S2 because tunneling Dirac fermions experience the direct bandgap near 
K- valley of M0S2, as aforementioned. 
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FIG. 4. Graphene/m-MoS2/graphene as a spin- filter for hole tunneling. Spin- 
dependent current density for a,b, electron and hole tunneling, respectively.through a graphene/m- 
MoS2/graphene heterostructure. Both panels exhibit spin-unpolarized currents (blue solid lines) in 
comparison with spin-dependent currents. Black and red solid lines indicate tunneling current den- 
sities for spin-up and down carriers, respectively. Green dashed lines represent the spin-polarization 
of tunneling current densities. 



Our results of spin-dependent tunneling currents through a magnetic M0S2 (m-MoS2) are 
shown in Fig. IHfor different kinds of carriers, electrons and holes. Here, the spin-polarization 
of tunneling current is calculated as Pj = {jup — jdown) / {jup + jdown)- In Fig- 111 one can see 
the differences in the spin-dependent features between electron and hole tunneling currents. 
While the hole tunneling current is almost perfectly spin-polarized, the electron tunneling 
current is weakly spin-polarized. This is due to the fact that the spin-splitting in the valence 
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band is about 50 times larger {5y ~ 150 meV) than the conduction band {5c — 3 mel^)|lO|. 
The small spin-splitting energy near K-valley in the conduction band leads to the weak 
spin-polarization of the electron tunneling current; Pj = 0.11 — 0.15. Meanwhile, due to the 
large spin-splitting near K-valley in the valence band, the hole tunneling current is almost 
perfectly spin-polarized, Pj = 0.99. Thus, the graphene/m-MoS2/graphene heterostructures 
provide a potential application in spintronic devices as a good spin filter. 

In summary, we have investigated the characteristics of the tunneling current through 
the graphene heterostructures with thin M0S2 layers. We have done calculations for vari- 
ous heterostructures; graphene/ Mo S 2 /graphene, graphene/MoS2/AGNR and graphene/m- 
MoS 2 /graphene. We show that the current on/off ratio of graphene /Mo S 2 /graphene is up 
to 10^ at room temperature. We also find that the graphene/MoS2 /graphene heterostruc- 
ture exhibits the carrier-dependent tunneling phenomena, which allows electrons to tunnel 
through the heterostructure, more easily than holes. Furthermore, we show that the current 
on/off ratio can be enhanced up to 10^ by replacing a two-dimensional graphene electrode 
with a GNR electrode. Finally, we propose a novel utility, the spin-polarized tunneling cur- 
rent, in the case that there is the splitting between different spin states in the thin M0S2 
layer. We show that the different spin-splitting energy between the conduction and valence 
bands of M0S2 makes the hole tunneling current almost perfectly spin-polarized whereas 
the electron tunneling current is partially spin-polarized. The graphene/m-MoS2 /graphene 
heterostructure may acts as a spin-filter for holes, which is a crucial building block of future 
spintronic devices. Our findings may not only offer an advance in research on the verti- 
cal graphene heterostructures with atomically thin insulating layers but also contribute to 
graphene-based spintronics. 
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